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Abstract 

The conditions under which relaxation dynamics in the presence of quenched-in disorder lead to 
rate-independent hysteresis are discussed. The calculation of average hysteresis branches is reduced 
to the solution of the level-crossing problem for the stochastic field describing quenched-in 
disorder. Closed analytical solutions are derived for the case where the disorder is characterized by 
Wiener- Levy statistics. This case is shown to be equivalent to the Preisach model and the associated 
Preisach distribution is explicitly derived, as a function of the parameters describing the original 
dynamic problem. 

Introduction 

There are many physical situations in which thermodynamic equilibrium is difficult to reach. 
The system may be trapped in long-living metastable states, and over ordinary time scales it may 
appear as permanently occupying a non-equilibrium state dependent on past history. Under these 
conditions, the relevant interpretative scheme is rate-independent hysteresis rather than equilibrium 
thermodynamics. In magnetism, micromagnetics and Landau-Lifshitz-Gilbert equations provide in 
principle the theoretical framework for the description of magnetization processes. Magnetic 
hysteresis should emerge from this description after appropriate averages over structural disorder. 
In its most general form, this problem is of discouraging difficulty. In this paper, we limit the 
discussion to a much simpler situation, that of a scalar system governed by relaxation dynamics of 
the form 

y^ = H(t)-H F (X) (1) 
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where y > is the friction constant, H(t) is the time-varying driving field and X(t) is the system 
response. The internal structure of the system is described by the free energy gradient H F (X). Due to 
structural disorder, H F (X) will fluctuate at random as a function of X. 

Neel [1] was the first to develop the concept of random energy landscape, in the framework of an 
explanation of the Rayleigh law on a statistical basis. Kronmuller and Reininger [2] started from 
this concept, investigating the properties of the derivative of the energy profile Hp(X) on the basis of 
the local fluctuations of the defect concentration, supposed to act as pinning centers. In our 
approach we will not enter in any detail concerning the physical origin of the fluctuations in the free 
energy gradient. 

We shall discuss the problem in the limit of low temperatures, where thermal agitation can be 
neglected and Hit) only is the relevant driving force. 

The hysteresis properties of the model emerge from the study of the trajectories H(X) of the 
system in the (X,H) plane. In particular, one finds that the rate-independent trajectories obtained in 
the limit \dH/dt\ — > are given by the sequence of points at which H F (X) reaches progressively 
higher field levels for the first time. On this basis, one can reduce the problem of the calculation of 
the average "magnetization curve" X(H) of the system to the solution of the level-crossing problem 
for the stochastic process Hp(X). By exploiting the mathematical results known in this field, we 
shall obtain an analytical solution for X(H) in the case where quenched-in disorder follows Wiener- 
Levy statistics. Remarkably, this analytical solution turns out to be equivalent to Preisach-type 
hysteresis. This result represents a limited but important example of derivation of macroscopic 
hysteresis properties from microscopic, non-equilibrium dynamic considerations. 

Hysteresis and the level-crossing problem 

The basic aspects of the dynamics described by Eq.(l) are made clear by considering the phase 
portrait of Eq.(l) in the (X,H) plane when the field increases at a given constant rate dH/dt (Fig. 1). 
Let us consider in particular the rate-independent limit dH/dt — > 0. In this limit, the trajectory 
obtained by integrating Eq. (1) is made of a sequence of stable states characterized by H = H F (X), 
connected by horizontal unstable branches (Barkhausen jumps) where dHldX = and dXJdt 
=(l/y)(H-H F (X)). This description can be recast in the following more general form, particularly 
useful when H F (X) is a stochastic process: the rate-independent trajectory under increasing field 
consists of the sequence of points H = H F (X) at which H F (X) reaches progressively higher field 
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levels H for the first time. This creates a direct connection between the trajectory and the statistics 
of the level-crossing problem for the stochastic process H F (X). 

Let us restrict the analysis to the case where H F (X) is a Markovian stochastic process. This 
means that the evolution of the process under given initial conditions, e.g., H F = Ho at X = Xo, 
depends on these conditions only and not on the details of the process for X < Xo. Let us consider 
the process in the interval (H_, H+ ), with H. < Ho < H+, as shown in Fig. 2. Level-crossing is 
studied by following the process in its diffusion away from H , until it reaches one of the levels H. 
or H+ for the first time. Here the process is stopped and removed. Mathematically, this feature is 
dealt with by introducing absorbing boundary conditions at H F = H. and H F = H+. The case relevant 
to the study of hysteresis is the one where a realization of the process starting at X = X. arbitrarily 
close to the lower boundary H. but not absorbed by it, crosses the upper boundary H+ at X = X+, 
with X+ > X. (see Fig. 2b). The crossing value X+ is a random variable, described by a certain 
conditional probability distribution, normalized to unity, which we denote by p LC (X + ,H+ I X.,H. ). 
There are methods discussed in the literature to derive this distribution [3, 4]. For the moment, we 
shall simply assume that it has been calculated and is known for a given process. 
Let us discuss how these results apply to the problem of hysteresis. First of all, notice that the 
system of Fig.l obeys the property known as return-point memory. This property can be defined as 
follows. Suppose H is decreased from large positive values down to Hi, then it is increased up to H 2 
and finally it is decreased again. If return-point memory holds, when H reaches again the value Hi 
the system returns back to the exact state it occupied when the field Hi was reached for the first 
time. It has been proven [5] that return-point memory holds for the individual trajectories of Fig. 1. 

Suppose we consider the trajectory of Fig. 1 when the field is decreased from some large 
initial positive value ( positive saturation ) down to H = Hi, where it is reversed and increased. Let 
X] be the value of X at the reversal point. The distribution p+(X,H ;Hi ) of the system response X 
after the reversal is controlled by the level-crossing probability p LC (X,H I Xi,Hi ) (see Fig. 2b, with 
H+ = H, H. = Hi) and by the distribution pi(Xi ; Hi) of the reversal value X], according to the 
expression 

X 

P+ {X,H;H l )=\p LC {X,H\X l ,H l )pXx i ;H l )dX l ; H > H, (2) 

Let us now suppose that the field is increased up to H = H2, X = X2, and then reversed again. Let 
P2(X 2 ;H 2 ,Hi) be the distribution of X at the new reversal point. p 2 (X 2 ;H2,Hi) coincides with p+(X 2 ,H 2 
;Hi) and can be calculated from Eq.(2). Each point (XH) of the new descending branch is a 
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potential upward reversal point. Therefore, because of return point memory, it is governed by an 
equation identical to Eq. (2), in the field interval (H,H 2 ). This means that the distribution p.(X,H 
\H 2 ,Hi) associated with the descending branch will be the solution of the inverse problem 

p 2 (X 2 ;H 2 ,H l ) = \p LC (X 2 ,H 2 \X,H)p.(X,H;H 2 ,H l )dX; H X <H<H 2 (3) 

This procedure can be continued to obtain p 3 (X 3 ;H 3 ,H2,H 1 ) = p.{X 3 ,H 3 ;H 2 ,Hi), and then 
p+{X,H;H 3 ,H 2 ,Hi) and so on. Any desired higher-order reversal branch can be generated and studied 
this way. The average system response along the generic n-th order branch will be 

X ± (H;H n ,...,H 2 M i ) = ]x P± (X,HM n ,...,H 2 ,H l )dX (4) 

where the "+" ("-") subscript denotes an ascending (descending) branch. We see that in general the 
higher-order branches will depend explicitly on the entire past sequence of reversal fields. 



Wiener-Levy statistics 

Let us apply the previous analysis to the particular case where the free energy gradient H F (X) 
has the form 

H F {X) = - + H P {X) (5) 

In Eq. (5), the term X/% represents some parabolic potential well ensuring the large-scale stability of 
the system, and H P (X), the pinning field, is a random process with zero mean, which we assume to 
coincide with the Wiener-Levy process. In particular, this means that <ldflpl> = 2AdX, with A > 0. 
Pinning fields with Wiener-Levy statistics are of definite physical interest. They were 
experimentally observed in ferromagnetic single crystals [6, 7], and their introduction leads to a 
satisfactory general interpretation of the Barkhausen effect [8]. 

The solution of the level-crossing problem for the Wiener-Levy process with drift (Eq. (5)) is 
known [4, 9]. One finds that Eq. (4) takes the simple analytical form 
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X ± {H;H n ,...,H 2 ,H ] )=(X n )±2A X : 



2A X 



2AX 



(6) 



A remarkable property of this solution is that it can be mapped onto the Preisach model of 
hysteresis. This derives from the fact that Eq.(6) obeys return point memory, and another property, 
known as congruency, which represent the necessary and sufficient conditions for the description of 
a given system by the Preisach model [9, 10]. The congruency property can be defined as follows. 
In a system obeying this property, reversal branches starting from identical reversal fields are 
geometrically congruent, whatever the past field history. Congruency does not hold in general for a 
system described by Eqs. (2)-(4). However, it holds in the particular case of Wiener-Levy statistics. 
In Eq.(6), in fact, the influence of past history is all lumped in the average value <X n > at the last 
reversal point. Therefore, all reversal branches have the same geometrical shape and merely differ 
by a shift, controlled by the values of H n and <X n >. This proves the congruency property. 

The conclusion is that there exists a Preisach distribution p(h c ,h u ) (with h c and h u representing local 
coercive fields and local interaction fields, respectively) that describes the average hysteresis 
properties of the dynamics in Wiener-Levy-type disorder. This distribution can be determined by 
first noting that, according to Eq. (6), any generic hysteresis branch depends on the difference H-H n 
only. It is known that, when this is the case, the Preisach distribution will depend on the local 
coercive field h c only. The dependence on h c is then obtained by calculating from Eq. (6) the 
derivative d 2 X + /dH 2 = (%/2) p[(H-H n )/2]. The result is [11]: 



p ( h 2 xctiix-1 X = A_ (7) 

c Ax sh 2 * ' Ax 



p(h c ) is normalized to unity, j p(h c )dh c = 1, and can be interpreted as a distribution of local pinning 
fields. Figure 3 shows a set of loops of increasing amplitude calculated from Eq.(6), or, 
equivalently, from standard Preisach integrals based on Eq.(7). The analytical derivation was tested 
by generating the same loops also through computer simulations. Individual random trajectories of 
the type shown in Fig. 1 were obtained from individual realizations of the process Hf(X) and were 
averaged over many subsequent realizations. The results of the simulations differed by less than 1% 
from the analytic results of Fig. 3. 
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Conclusions 

We have seen that, whatever the statistics of the pinning field of Eq. (5), return-point memory 
will hold under all circumstances. Congruency, on the other hand, was obtained only in the 
particular case of Wiener-Levy statistics, and in general will fail for other types of random fields. 
This raises the interesting question of what type of property could play the role of congruency under 
more general conditions. Answering this question would show the direction in which one should 
generalize the Preisach model in order to approach a universal description of hysteresis and would 
give a natural way to introduce new classes of hysteresis models. 
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Figure Captions 

Fig. 1 - Free energy gradient H F (X) and corresponding individual trajectory H(X) in the rate- 
independent limit \dH/dt\—>0. The state (Xj,Hi) is an upward reversal point. 

Fig. 2 - Schematic representation of the level-crossing problem in the field interval (H.,H+), starting 
from a given initial condition (Xo,Ho) (a), and in the limit Ho—> H. (b). 

Fig. 3 - Hysteresis loops calculated from Eq.(6). In the inset is shown the Preisach distribution p(h c ) 
(Eq.(7)). 
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